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V578 Mon is an eclipsing binary system in which both stars have masses above 10 Mq de- 
termined with an accuracy better than 3%. It is one of only five such massive eclipsing bina- 
ries known that also possess eccentric orbits and measured apsidal motions, thus making it an 
important benchmark for theoretical stellar evolution models. However, recently reported de- 
terminations of the radii of V578 Mon differ significantly from previously reported values. We 
reanalyze the published data for V578 Mon and trace the discrepancy to the use of an incorrect 
formulation for the stellar potentials in the most recent analysis. Here we report corrected radii 
for this important benchmark eclipsing binary. 



Introduction 



In a recent major review of eclipsing binary (EB) stars, iTorres et al.l ( 20101 ) compiled a sample of 95 
"benchmark grade" EBs having masses and radii measured with an accuracy of better than ±3%. This 
sample of EBs is of funda mental importance t o a broad range of applications in stellar astrophysics. For 
example, using this sample, [Torres et al.l (|2010[ ) derive empirical relations that permit the determination of 
fundamental properties, such as mass and radius, for any star from its directly observable properties, such as 
effective temperature, surface gravity and metallicity. This sample of EBs is moreover important for testing 
the basic predictions of theoretical stellar evolution models, which are in turn used across many areas of 
research, including star formation, stellar evolution, and exoplanets (because the inferred properties of the 
planets generally depend on the assumed properties of the host stars) . 

One of the 95 EBs included in the lTorres et al.l ( 20101 ) compilation is V578 Mon. This EB is remarkable 
by virtue of being one of only nine EBs in the compilation where both stars have masses greater than 10 
M f7), and of these, is one of only five to have an ecce ntric orbit and therefore a measurable apsidal motion 



( Hensberge et al.ll2000t Pavlovski fc Hensbergd 12005 ). Indeed, the apsidal motion of V578 Mon was only 
recently measured (JGarcia et al.ll2011[ ). 



We have identified an error in the calculation of the stellar radii reported in the iTorres et al 



(201 



3) 



compilation. In this brief report, we trace the source of this error to the use of a formulation for the stellar 
potentials that does not correctly include the system's orbital eccentricity. We then provide updated radii 
calculated from the correct formulation for the potentials. 



2. Calculation of Stellar Radii for V578 Mon 



Torres et al.l (|2010l) reported the radii of V578 Mon to be Rp =5.149 ± 0.091 and Rs =4.21 ± 0.10, 
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where Rp is the radius of the more massive primary star and Rs is the radius of the less massive secondary 
star. These rad i i diffe r from the previously reported values of Rp =5.23 ± 0.06 and Rs =4.32 ± 0.07 by 
Hensberge et al.l ( 20001 ). While these differences are not highly statistically significant given the uncertainties 



quoted in the two studies, the difference is slighly more than la and any differe nce is unexpected sinc e 
Torres et al.l ( 20101) recalculated the radii from the same data originally presented bv lHensberge et al.l (2000). 



As is standard for the determination of EB radii, the radii are calculated from the volumetric radius 
which is defined as R = {-^VY^^ ^ where V is the volume of the star. Therefore we suspected the difference 
may lie in different approaches to calculating the stellar volumes, which are themselves derived from the 
ste llar potential surfaces. Since [Torres et al.l (|2010[ ) adopted the same stellar potentials originally reported 
bv lHensberge et al.l ( 20001 ). the different radii reported by the two authors is likely to be in the formulation 
used to transform the potentials to volumes. 

Torres et al. ( 20101) com putes V using the stellar potentials, 17, and the stellar mass ratio, q = M2/M1, 
from lHensberge et al.l ( 20001) assuming a circular synchronous orbit. More formally, V = V{^, q) if the orbit 
is circular (e = 0), and if the stars rotate synchronously with the orbit ( F^ = Fy = 1. 0,where F is the ratio 



of the stellar angular velocity to the orbital angular velocity). Following iKopall (Il989l ): 
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where v = cos0, A = cos (^ sin 0, 9 and are instantaneous direction angles between the two stars, and r 
is the distance to the surface of the star from the star's center in the (6*, 0) direction and in units of the 
instantaneous separation between the two stars. Through several complex mathematical operati ons, equatio n 
|l]can be used to uniquely define V = V{Vl, q) and therefore R = -R(il, q) (see Equation 2.18 in lKopallll989l) . 



However, V578 Mon has an ecce ntric (i.e., non-circular ) orbit with pseudo-synchronous rotation, where 
e = 0.0867, Fi = 1.13 and F2 — 1.11 ( Hensberge et al.ll2000 ). Therefore the R cannot be uniquely determined 
as a fu nction of fl a nd q alone. Instead of using equation [U it is necessary to use VI defined for non-circular 
orbits (|Wilsonlll97"9l ): 
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where b = D j a is the instantaneous separation between the two stars normalized to t he semi-major axis . 
Equation H thus defines V = V(Sl,q,F,b) and t herefore R = R(n q,F 6). Whereas l^rres et al.l |2010l ) 



calculated the V5 78 Mon radii via eq uation [H iHensberge et al.l (|2000h uses the Wilson-Devinney code 



m 



(jWilson fc Devinnevlll97ll : IWilsonlll979[ ) which correctly incorporates equation |2| in computing the radii. 



As an independent check on the Wilson-Devinney based calculation of the radii, we compute Rp and 
Rs by solving equation |2| nume r ically for r{9,(j)) in unit steps of 9 and (f>. We set ili, fl2, Fi, F2, q and e 



to values from lHensberge et al.l (2000|). For noncircular orbits, H. = n{t) since the instantaneous separation 



between the two stars varies over the orbit. We compute th e volumetric radius at p eriastron, setting the 
instantane ous sepa.ration S = 1 — e, since the ft reported bv IHensberge et al.l (|2000[ ) are the potentials at 
periastron ( Wilsonll 19791) . We compute V from r{9,(j)) discretely as: 



where A(f> and A9 arc discrete steps in (f> 



3 

and 9. 



(3) 



From equation 13 w e compute Rp = 5.24 and Rs = 4.33, which are essentially identical to the radii from 
Hensberge et al.l (|2000l ). Finally, in the same fashion we also solve for r{9, (f>) for the case of a circular orbit 



(equation [T]), and compute V using equati on [3l In this c ase, w e find Rp 
essentially identical to the radii reported in [Torres et al.l ( 2010r ). 



5.15 and -Re = 4.21, which are 



Thus, it is clear th at the differences in the radii re ported for V578 Mon between iHensberge et al.l ( 20001 ) 
and [Torres et al.l (|2010l ) arises from an incorrect use by [Torres et al.l (|2010( ) of the formulation for the stellar 
volumes assuming circular, synchronous orbits. 



3. Summary 



We have shown that the radii report ed for V578 Mon by Torres et al.l ( 20101) are incorrect, with an error 
of ~2% relative to the expected values (IHensberge et al.ll2000l). W hile small, the error is unexpected since 
Torres et al.l (|2010l ) used the same data as IHensberge et al.l ([20001), and moreover the highest quality radius 
determinations of EBs should achieve ~1% precision. We have traced the error to the use by those authors 
of the formulation for the stellar potential surfaces that is correct only in the case of circular, synchronous 
orbits. However, V578 Mon possesses an eccentric, pseudo-synchronous orbit, and thus its volumet ric radii 
must be calculated from the formulation for the stellar potentials that correctly accounts for this ([Wilson 



19791) 



Radii for V578 Mon were originally reported bv IHensberge et al.l (2000), who use the Wilson-Devinney 
code which incorporates the correct formulation for t he stellar potent i als in an eccentric orbit (equation 
[2]). We have performed an independent check on the lHensberg;e et al.l (|2000l ) radii by solvin g equation 
numer ically: we recover radii identical to IHensberge et al.l ( 20001 ). We are able to recover the iTorres et al.l 



i 



(|2010[ ) radii when we solve the equations assuming a circular, synchronous orbit. 



We have checked that the use of the erroneous formulation for the potential in lTorres et al.l (|2010l ) does 
not significantly alter the radii for any other systems reported in that paper; fortunately only V578 Mon is 
sufficiently massive and eccentric to have been affected. A forthcoming paper (Garcia et al., in prep.) will 
report improved radii of V578 Mon ta king into account an updated eccentricity of e = 0.07755 and apsidal 
motion to = 0.07089 deg cy cle 



(|2011l ) . In the meantime we recommend the use of the 



radii from IHensberge et al.l (|2000l ) which we have corroborated here. 
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Tabic 1. V578 Moii Radii Comparisons 







Rp 


Rs 


Hcnsbcrge et al. f2000~) 


5.23 ±0.06 


4.32 ±0.07 


Torres ct al. C2010) 


5.149 ±0.091 


4.21 ±0.10 


This Study 


5.24 


4.33 





